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We present numerical results in SU{2) lattice gauge theory for the space-space and 
time-time components of the gluon propagator at equal time in the minimal Coulomb 
gauge. It is found that the equal-time would-be physical 3-dimensionally transverse gluon 
propagator D^^{k) vanishes at /c = when extrapolated to infinite lattice volume, whereas 
the instantaneous color-Coulomb potential D44(A;) is strongly enhanced at /c = 0. This has 
a natural interpretation in a confinement scenario in which the would-be physical gluons 
leave the physical spectrum while the long-range Coulomb force confines color. Gribov's 
formula D^^{k) = (|fc|/2)[(P)^ + M^]^/^ provides an excellent fit to our data for the 3- 
dimensionally transverse equal-time gluon propagator D^^{k) for relevant values of k. 
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1. Introduction 

Wilson's lattice gauge theory provides a regularized formulation of gauge theory that is 
manifestly gauge invariant, and numerical simulations do not require gauge fixing. However 
gauge fixing on the lattice is advantageous to gain control of the critical or continuum limit, 
for this makes available the strong results of gauge- fixed continuum renormalization theory. 
For example one may prove in continuum renormalization theory that a certain quantity, 
such as the running coupling constant, defined in terms of gluon correlation functions in 
some gauge, is finite when the cut-off is removed. Then the corresponding lattice quantity, 
defined in the corresponding lattice gauge, should be finite in the critical limit, and it 
becomes of interest to make a numerical determination of that quantity. Moreover one 
may determine by numerical fit the location of the poles of propagators which, according 
to the Nielsen identities, is independent of the gauge parameters [1] . Finally, if one has in 
hand a confinement scenario in a particular gauge, then it is possible to test its predictions 
numerically for gauge-fixed quantities. Although the scenario may look quite different in 
different gauges, nevertheless any one of them provides a valid perspective. 

Previous numerical studies of the Coulomb gauge were reported in [2] and [3]. We 
present here a numerical study of the gluon propagator in SU{2) lattice gauge the- 
ory, without quarks, in the minimal Coulomb gauge (defined in the Appendix). Sim- 
ulations have been done at /? = 2.2 for 9 different lattice volumes V = L^, with 
L = 14, 16, 18, 20, 22, 24, 26, 28 and 30. (A total of 2420 configurations have been 
generated, from 50 configurations for 30^ up to 600 for 14^.) The procedure is to first 
equilibrate ungauge-fixed configurations U according to the Wilson action using a hybrid 
over-relaxed algorithm [4] . Then statistically independent configurations U are gauge fixed 
to the minimal lattice Coulomb gauge by a minimization that is effected using a stochas- 
tic over-relaxation algorithm described in [5] with accuracy {{diAi)'^) < 10~^^, where the 
average is taken on each time slice separately and diAi is defined in Eq. (A. 4). Finally 
the components D^^{k) and D44^{k) of the equal-time gluon correlator are evaluated. The 
lattice Coulomb gauge is more easily accessible to numerical study than the Landau gauge 
because each time-slice contributes separately to the numerical average which, for a lattice 
of volume 30"^, gives a factor of 30 gain. The total computer time devoted to this project 
so far is about 500 days on a 500 MHz ALPHA work-station^. 

^ We thank Jorge L. deLyra for kindly providing us with access to the cluster of ALPHA 
work-stations at the Department of Mathematical Physics (DFMA) of the University of Sao Paulo 
(USP). 



1 



The results provide a test of the confinement scenario that was originally proposed 
by Gribov [6] and elaborated in [7]. The confinement scenario is particularly transparent 
in the Coulomb gauge because it is a physical gauge in the sense that the constraints are 
solved exactly, including Gauss's law DiEi = diEi + [Ai, Ei] = pqu, and the Hilbert space 
has positive metric. Here E^ is the color-electric field, Di is the gauge-covariant derivative, 
and pqu is the color-charge density of quarks. Moreover the gauge fixing, described in the 
Appendix, is done independently within each time-slice so that the equal-time Euclidean 
and Minkowskian of the correlation functions are identical. 

In the Coulomb gauge, the 3- vector potential Ai is transverse, diAi — 0, so A^ — Af^. 
Gauss's law is solved by 

E, ^ El' - dicl>, (1.1) 
where the color-Coulomb field is given by 

= M-Vcoui- (1.2) 

Here pcoui = Pqu — [^f' color-charge density of the dynamical degrees of freedom, 

and M = M{A^'^) = —Di{A^')di is the 3-dimensional Faddeev- Popov operator. 
In this gauge, the Hamiltonian (without quarks) is given by 

H = (2^2)-l J £^ (£;2 ^ ^2) 

= {2gl)-^ J d^x {E'^^ + B^) + {2gl)-^ J d^x d^y p,oui{x) V{x,y;A'') peoui(2/) • 

(1.3) 

Here 

V{x,y;A"^) = [M{A'T\-d^)M{A'T%,y (1-4) 

is a color-Coulomb potential-energy functional, depending on A*^, that acts instantaneously 
and couples universally to color charge. The continuum form of the Coulomb Hamiltonian, 
with proper attention to operator ordering is given in [8] and the lattice form in [9] . For the 
minimal Coulomb gauge, this Hamiltonian is supplemented by the boundary condition that 
the wave functionals \1/(A*'') are restricted to the Gribov region. (This is explained below 
and in the Appendix.) It was proposed in continuum theory in [6] and in lattice theory in 
[10], that this restriction may be imposed by use of an effective action or Hamiltonian 

H,s = H + {2gl)-^ j d'x Af{-VY'Af . (1.5) 



As a result, the energy of a gluon of momentum k gets modified to E'^ik) = k'^ + {k'^)~^M'^, 
as one sees from the quadratic part of ffeff • One obtains for the equal-time 3-dimensionally 
transverse would-be physical gluon propagator the approximate expression 

D''(k) = (27r)-i / dk4 ^ = -^-^ (1.6) 



This quantity is defined by 



D''(k) = ^ . (1.7) 

2[(A;2)2 + M4]V2 



Dfj{k) = {Sij-k,kj)D''{k) , (1.8) 



where Dlj{k) is the Fourier transform of 



DlUx-y) = {Af{x,t)Af{y,t)) . (1.9) 



In the absence of an estimate of corrections, one does not know how accurate (1.7) may 
be. However it was proven [11] that the lattice gluon propagator D^^{k) at infinite spatial 
lattice volume must indeed vanish at ^ = 0, 

limD^'^(^) = 0, (1.10) 

although the rate of approach of D*'^(0, L) to 0, as a function of L, was not established, nor 
was it determined whether the renormalized gluon propagator also shares this property. 
The accuracy of (1.7), and the crucial question of the extrapolation to large L of D^^{k, L) 
are addressed in the numerical study reported here. 

We have also evaluated the 4-4 component of the gluon propagator 

i?44(x, t) = {A^ix, t)A40, 0)) (1.11) 

at equal-time. In the minimal Coulomb gauge D4^4^{x,t) is given by [12] 

DM(x,t) = V(x)S{t) + P{x,t) , (1.12) 

where t = X4 is the Euclidean "time" and P{x,t) is a non-instantaneous vacuum- 
polarization term. This gives 

J dt D44{x, t) = V{x) + o(e), (1.13) 
3 



where o(e) vanishes with e. We call V{x) the color-Coulomb potential. In momentum 
space (1.12) reads 

D44{k, k4) = Vik) + P{k, /C4) , (1.14) 
where lim/,^^00 P{k, k^) = 0. In dimension d < 4, V{x) coincides with 

Voix-y) = {V{x,y;A'^)). (1.15) 

However in d — 4 dimensions there is a mixing of V{k) and P{k, /c4) associated with 
divergences, and V{k) differs from Vo(fc) by terms of the form CngQ^/k^ in each order of 
perturbation theory, as explained in detail in [12]. 

Stated simply, the confinement of color is caused by the predominantly long range of 
the color-Coulomb potential V{x), corresponding to an enhancement of V{k) at low |^|. 
Note however that V{x) is not the gauge-invariant energy eigenvalue of the quantum state 
of infinitely massive separated quarks. Nevertheless it is an important quantity. It may 
be used as an order parameter for color confinement [7], and it is the starting point for 
calculations of the ground-state wave-functional [13], [14] and [15]. 

The rather surprising and counter-intuitive vanishing of D^^{k) at k = 0, and the 
enhancement of y(/c) at /c = in the minimal Coulomb gauge are both caused by the 
Gribov horizon. This is a boundary in the space of configurations Af(x) defined by the 
condition that the Faddeev-Popov operator be positive, M{A^^) > 0. The Gribov horizon 
represents the points where the lowest eigenvalue Xo{A^^) of M(A*'') first goes negative. 
As shown in the Appendix, all configurations that contribute to the Euclidean functional 
integral in the minimal Coulomb gauge are constrained to lie within the Gribov horizon. 
Because of entropy considerations (see Refs. [6] and [10]), the Euclidean probability gets 
concentrated near the horizon where the color-Coulomb interaction energy V(^*'^), Eq. 
(1.4), diverges. This causes an enhancement of the instantaneous color-Coulomb poten- 
tial V{x). At the same time it places very severe bounds on the magnitude of the low 
momentum components of the gluon field. ^ 

The Coulomb gauge does not offer any particular advantage for perturbative calcu- 
lations and renormalization. However it is the finite limit of renormalizable gauges [19]. 

^ A similar confinement scenario for the Landau gauge was proposed in [16] and [17]. It has 
also been verified numerically in the Landau gauge that typical (thermalized and gauge-fixed) 
configurations lie very close to the Gribov horizon [18]. 
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A valuable feature of this gauge is that the time-time component of the gluon propaga- 
tor D44(k,k4) is independent of both the cut-off A and the renormalization mass fj, [7]. 
This holds separately for its instantaneous part V{k). Thus the minimal Coulomb gauge 
allows us to introduce a running coupling constant by ^coui(l^l) ~ const k'^V{k). The 
proportionality constant is determined by the condition that in the large-momentum or 
weak-coupling regime g^ouiil^l) satisfy the standard renormalization-group equation, 

1^1 %^ = /3coul(^coul) = -(Mcoul + ^l^coul + • • •) (1.16) 

o\k\ 

where, for SU{N) gauge theory without quarks, bg = (47r)-2lliV/3, bi = {47r)-'^34N'^ /3. 
The proportionality constant is calculated in [12], with the result that for SU{N) gauge 
theory without quarks 

k'V{k) = g ^e'ould^l/Acoul) , (1-17) 

and more generally, with Nf quark flavors 

^'^^^^ = nN^-2Nf ^couid^l/Acoui) . (1.18) 

Here Acoui cx: Aqcd is a finite QCD mass scale, characteristic of the Coulomb gauge, such 
that asymptotically in the weak-coupling regime 

P = Al,, exp[(6o^e'oul)"'] (^'O^c'oul)'^ , (1-19) 

where r = bi/bl. The ratio Acoui/Aqcd may be obtained from a 2-loop calculation [12]. 
We conclude that in the minimal Coulomb gauge the running coupling constant of QCD 
may be obtained from a numerical determination of the equal-time 2-point function D44, 
whereas in other gauges it must be obtained from a 3-point function. The running coupling 
constant ^fcoui that we have introduced is the QCD analog of the invariant charge in QED 
that is defined in terms of the transverse part of the photon propagator in a Lorentz- 
co variant gauge. 

2. Results 

We evaluate the space-space gluon propagator 

L 3 3 

w 



t=i fi=i b=i 



*="■="-' (2.1) 

L 3 3 ^ ^ 



t=l ij.= l 6=1 
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and the time-time gluon propagator 



where 



t=l 6=1 



(2.2) 



+ 



E ^'^t, sin {9 



X 



(2.3) 



and the lattice gluon field A'l^^ is defined in Eq. (A. 5). Here, as usual, we have defined 
ki = 2sin(0i/2), for — tt < 9i = 27rni/L < tt, and integer n^. An average over the L time 
slices is included. In our simulations we consider only 3-momenta aligned along major axes 
Oi = (0, 0, 27rn/L). Notice that D*='(0) is not given by D^'{k) at k = 6. The difference is 
due to the Coulomb gauge condition — the continuum-like condition, Eq. (A. 4) — which 
in momentum space reads 

E^.^.-t. = 0, (2.4) 

i=l 

■ — ' — * 

where A^^^ is the three-dimensional Fourier transform of the gauge field A^u. \i k ^ 
(0, 0, 0) only two of the three Lorentz components of A^^^ — and therefore of A^ti — are 
independent. This explains the factor 6 (instead of 9) in the definition of D^^ik). 

Fig. 1 shows D*^(A;) and Z>44(A;) as a function of k on the same logarithmic plot for 
the lattice sides L = 28 and 30. The qualitative behavior is quite different in the two cases. 
Whereas -D44(fc) grows strongly at low fc, by contrast D*^(fc) turns over and decreases at 
low k. 

(a) Analysis of D*''. Figs. 2, 3 and 4 show our data for D^^{k) for the nine lattice 
volumes considered. To parametrize these data, we were guided by the Nielsen identities 
[1]. They tell us that the poles of the gluon propagator D^^{k, k^) are independent of the 
gauge parameters.^ In particular, for the class of gauges defined by A ^4^4 -|- diAi = 0, 



^ Strictly speaking, the Nielsen identities have been established only for Faddeev-Popov type 
gauge fixing. This does not include the minimal Coulomb gauge because the gauge fixing is done 
by a minimization procedure that is not describable by a local 4-dimensional action. However the 
minimal Coulomb gauge may be obtained as a limiting case of a local 5-dimensional quantum field 
theory that describes stochastic quantization with stochastic gauge fixing, as has been discussed 
recently in [16] and [17]. The Nielsen identities may be extended to the 5-dimensional formulation. 
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which interpolate between the Landau gauge, A = 1, and the Coulomb gauge, A = 0, the 
poles are independent of A. In the Landau gauge, the poles occur at = k'^ + kl = — m^, 
by Lorentz (Euclidean) invariance, and thus also in all these gauges, by virtue of the Nielsen 
identities. We have made a 2-pole fit with poles at mf and m|, which may be either a 
pair of real numbers or a complex conjugate pair.^ According to [11], the propagator in 
the minimal Coulomb gauge vanishes at k — 0. We use this condition to fix the residues 
to within an over-all normalization 



D'\k,k^)^C (P)"-^ 



1 / kl + ml kl + ml 



k'^ + kl + ml k^ + kl + ml'' 

-C (P)'^ ( ^ ^ ] 

^ ' ^P + kj + ml P + kl + my ' 



(2.5) 



where ct is a fitting parameter. This formula does not reproduce the correct asymptotic 
behavior at large momenta ~ (/c^ + ^1)"^ times logarithmic corrections. However, this is 
not a problem because our largest value of |^| is less than 2 GeV, and we are probably far 
from the ultraviolet regime.^ We wish to emphasize that our pole parameters — to the 
extent that they are a valid fit — are gauge-independent quantities that characterize the 
gluon. 

The equal-time part of this propagator is given by D^^{k) = {2tv)~^ J dO^ D^^{k, /C4), 
where k4 — 2 sin (6*4/2). This gives 

D'^ik) = -C' {Pr (A-'/^ - A^'^^) 

^ (P^a 4 + fei + /i2 (2.6) 

^ ^ Al/'Al/'{Al/' + Al/') ' 

where hi = P + mf and Ai = Ahi + hf for z = 1, 2. If we take a pair of complex conjugate 
poles, mf = X + iy and m^ = x — iy, we obtain 

(u2+ V^;2) 1/2 [(^^2 +4^2^2)1/2 1/2 ' [^■') 



^ According to the general principles of quantum field theory, the propagator of physical parti- 
cles should have poles only at real positive . However in the confined phase the gluon propagator 
may have singularities that correspond to unphysical excitations. 

^ In the Landau gauge [20] the gluon propagator reaches three-loop asymptotic scaling at 
momenta of about 5-6 GeV. 
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where u = 4(^^ + x) + (P + x)^ — and v = + + x). The case of a pair of real poles 
is obtained from this formula by taking a negative value of y^. For the case of a lattice of 
finite volume F = L^, we modify this formula to 

D^'{k) = Z [(P)'^ + r] (^^2 + 4^2^2)1/2 [(^2 + 4^2^2)1/2 + ^]l/2 " (2-8) 

The factor v[(w2 + 42/2^;2) 1/2^^] -1/2 ^ ^A^- A2)[^^'^iy{A\''^ + A^^)]-^ is slowly varying 

— * 

over the relevant range of k and parameter values, so a fit to (2.8) is a test of the simpler 
formula 

D-(k)=^' [(gr + '-l („.^4^.„.)./. - (2.9) 

Stated differently, (2.8) and (2.9) have the same singularities that are nearest to the origin. 
In the continuum limit we have u A(hP'-\-x) and v ^ 2^ and (2.9) is a lattice discretization 

of Gribov's approximate formula (1.7) provided that the fitting parameters have the values 
r = 0, a = 0.5, x = 0, with the identification y"^ = M^. It is intended report on the fit to 
(2.9) elsewhere, but preliminary indications are that it is comparable in quality to the fit 
to (2.8). 

For each lattice side L we have made a fit of the parameters z{L), r{L), a{L), x{L) 
and y^{L). By using Table 3 of [21] and by setting the physical string tension equal to 
^/a = 0.44 GeV we obtain that, for /? = 2.2, the inverse lattice spacing is = 0.938 GeV. 
This gives a = 0.21 fm, so that the largest lattice volume considered here, i.e. V = 30*, 
corresponds to (6.3 fm)*, the smallest non-zero momentum that can be considered for that 
lattice is equal to 0.196 GeV, while the maximum momentum value (for each lattice side 
L) is 1.876 GeV. The results^ are exhibited in Table 1, and the curves are plotted in Figs. 
2, 3 and 4. There is no a priori reason why a 2-pole fit should be accurate over the whole 
range of momenta considered. However the fit is excellent for all momenta /c^ and for each 
L. We have also checked that results in agreement with those reported in Table 1 are 
obtained if one considers only the data corresponding to /c^ < 2. 



^ The fits have been done using gnuplot; the errors represent 68.3% confidence interval. 
Similar results have been obtained using a conjugate-gradient method with errors estimated by 
a jack-knife method. Let us notice that we did not consider the correlations between different 
momenta when fitting the data. However, at least in the Landau gauge, the covariance matrix for 
the gluon propagator in momentum space is essentially diagonal [22] , and the value of the "naive" 
X^/d.o.f. is usually compatible with the value obtained using the full covariance matrix [23]. 
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We have extrapolated to infinite L the fitting parameters z{L)^ ol{L)i x{L), y'^{L) and 
the product {rz){L)\ in all cases we tried three different types of fitting functions, namely 
a+b/L^, exp(a)/L'' and a+log(l+6/L), and chosen the fit with smallest x^/d.o.f. Results^ 
are shown in Table 2 (first row) and plotted in Figs. 5-9. 

The reader will have noticed that the product {rz){L) extrapolates to (see Fig. 
9). This corresponds to a vanishing of D^^{0) at infinite lattice volume, in accordance 
with [11]. To check on this important point we have also fitted D*^(0,L) using the three 
fitting functions considered above. A good fit is provided by D^^{0,L) = exp(a)/L'', with 
a = 2.43(4), b = 0.50(1), xVd.o.f. = 0.69 and goodness-of-fit Q = 65.6% (see Fig. 10). ^ 
An indication of how reliable these fits are is the comparison of the power b in the fit of 
the product {rz){L) and of D^'^{0,L) which arc b = 1.1(1) and b — 0.50(1) respectively. 
We have also made a similar fit for r(L) (not plotted) and obtained b = 1.8(1). 

We next consider the k dependence of D^'^{k) at low k. For the power dependence 
parametrized by (P)'*, observe that a{L) extrapolates to a{oo) = 0.48(5), see Fig. 6. 
This agrees with Gribov's approximate formula (1.7), which gives a = 0.5. Moreover this 
value is consistent with the other rows of Table 2, by the method described below, so 
this result appears quite stable. Particularly striking is that, with a; = imposed, one 
obtains a = 0.49(1) and a = 0.51(1) respectively from the fourth and fifth rows of Table 2 
(explained below). 

Another striking feature of the fit is that y'^{L) is positive for all 9 values of L (see 
Table 1), corresponding to a pair of compiea; conjugate poles rather than a pair of real poles. 
Moreover, in all (L) is quite small compared to y{L). For example a;(30) = 0.06(6) 

and y{SO) = 0.88(5). The extrapolation to infinite L also strongly indicates a positive 
value y^(oo) > and a small, possibly zero, value for x{oo). Thus our data arc compatible 
with and perhaps suggestive of poles at purely imaginary — ± iy, in agreement with 
Gribov's formula (1.7) with y"^ = M^. 

Finally, we have fit the data for L = 28 and L = 30 using Eq. (2.8) with r = and 
with a low- momentum cut /c^j^^, namely considering only a range of momenta in which 
finite-size effects are negligible. Results are reported in the second and third rows of Table 

^ These fits have also been done using gnuplot. For the fitting function exp(a)/L'' the results 
have been checked with the exact minimizing formula (notice that this fitting function is linear in 
the coefficients a and h after taking the logarithm). 

^ When using the fitting function a + b/L'^ we obtain, both for {rz){L) and D*'^(0, L), a value 
of a that is zero within errors but a worse x^/d.o.f. 
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2 for two different values of ^min- Similar results are also obtained when the fixed value 
a; = is imposed (see the last two rows of Table 2 and Fig. 11). The values for a and 
obtained in this way are in good agreement with the values obtained by extrapolating 
a{L) and y'^{L) to infinite L (compare the first row of Table 2 with the other four rows of 
the same table). 

(b) Analysis of -D44. Eq. (1.13) shows that for lattice quantities we may make the 
identification V{x) — L>44(if), where D44{x) is the equal-time propagator, and similarly 
for their 3-momentum transforms, V{k) = D44{k). This allows us to use (2.2) and (1.17) 
to also define the lattice quantity fifcouil^) 

^9coui{k) ^ k^Vik) ^ PD,,{k) . (2.10) 

Fig. 12 shows our data for the running coupling constant fifcouil^)- continuum 
limit, its behavior at large momentum is governed by the perturbative renormalization 
group (1.16) and (1.19). For the fitting formula we modify (1.19) to 

= ALui eMibglui)-'] [{bgluiy + <bglj^]-' , (2.11) 

which implicitly defines fifcouil^^)- Here r = 102/121, and A^g^j, 6, 2;, a are fitting parame- 
ters whose significance we now explain. Naturally the parameter ■A^^.^^j sets the mass scale. 
For small g'^^^i, which corresponds to large k'^, this formula is dominated by the first term 
in the denominator, whereas for large (/^ouIj which corresponds to small fc^, it is dominated 
by the second term in the denominator. For small fi'^oui formula approaches (1.19), pro- 
vided that b = bo = ~ 0.046. However we are quite far from the continuum limit for 
i:>44(^), as indicated by the small value ((l/2)tr(t/4)) = 0.221(6) (for lattice volume 14"^),^° 
and we expect significant /9-dependence in the extrapolation to the continuum limit. More- 
over, for fixed /3, it has been found that different lattice discretizations of the gluon field 
lead to identical gluon propagators to within numerical accuracy, apart from the overall 
normalization [24]. We allow for this by taking the overall normalization of g^^y^x to be a 
fitting parameter. This requires putting an arbitrary normalization coeflftcient everywhere 
in front of g^^^x in (2.11), which is equivalent to replacing the fixed number 60 by the fitting 
parameter h. Of course, an extrapolation in /? to the continuum limit should give b — bo. 

1° Before minimizing Fver we find ((l/2)tr(?74)) = -0.0005(57). 
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For large g^^^i and small fc^, this formula approaches 

'>9L.= {^) ■ (2.12) 

Thus the parameter z sets the overall normalization in the strong-coupling or infrared 
regime, and a governs the strength of the singularity of ^^0^1 in the infrared limit. If 
the color-Coulomb potential V{k) is governed by a string tension at large distances then 
5'coui(^) ~ const/^^, which corresponds to a = 1. 

On a finite periodic lattice with finite lattice spacing and with fi'coui(^) defined by 
Tf5'coui(^) = k'^V{k), necessarily glQ^\(k) vanishes at = 0. Actually fi'coui(^) S®^^ ^ 
maximum value at fc^ fa 0.2 and goes to zero in the infrared limit (see Fig. 12). This 
unphysical behavior is a lattice artifact that also appears in studies of the running coupling 
constant using the three-gluon vertex [25]. In order to fit our data we have made a 
low- momentum cut at fc^jn = 0.5. The values of the parameters which we obtain^ ^ are 
A2^^j = 1.0(2), z = 1.2(1), h' = 0.18(2), b = 0.20(2), a = 1.9(3), with x^/d.o.f. = 0.44 and 
goodness-of-fit Q — 98.3% using the data for L = 28 and L = 30, where b' = ^b. We have 
checked that similar results are obtained with fc^j^ = 0.3 and k^i^ = 1.0 (see Table 3) but 
the resulting x^/d.o.f. is smallest for k^^^ = 0.5. 

The value of a ~ 2 corresponds to g'^^^i ~ const/ and V{k) ~ const/ at low 
momentum. The volume dependence of the data, and therefore of our fit, is quite weak, 
but one notices in Table 3 that the value of a decreases as the low momentum cut-off k'^^^ 
increases. This corresponds to an increase in the strength of the singularity of fifcouil^) 

— * 

A; = as finite-volume effects are reduced. However, as explained above, we must make 
an extrapolation in /? in order to arrive at any precise conclusion about the strength of 
the singularity in the continuum limit. Nevertheless, our data at finite P and L clearly 
indicate a color-Coulomb potential that is more singular than V{k) ~ const/ |/c| at low k. 

3. Conclusions 

We have used simple formulas to fit the data for the equal-time gluon correlators 
D^^{k) and 1^44 (/c) in SU{2) lattice gauge theory in the minimal lattice Coulomb gauge at 
j3 — 2.2. Our fits have the following features: 

The fit has been done using a conjugate gradient method with a numerical inversion of Eq. 
(2.11). Errors are estimated using a jack-knife method. 
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1. The equal-time would-be physical gluon propagator D^^{k,L) at zero momentum 
k extrapolates to in the limit of infinite lattice volume, i.e. D*^(0, oo) = 0. The rate 
of approach for lattice volume V = L"^ was fit by D^^{0, L) = CL~^, where h = 0.50(1). 
The vanishing of the gluon propagator D^^{k) at = is highly counter- intuitive, and the 
only explanation for it is the suppression of the low momentum components of the gluon 
field that is a particular feature of the Gribov horizon, that constitutes the boundary of 
configuration space. -"^^ 

2. Asymptotically at low momentum our fitting formula behaves like D^^(k) oc (A;^)°, 
with the value extrapolated to infinite lattice volume a — 0.48(5) (see Table 2). With the 
value a; = imposed and a low-momentum cut-off one obtains (from L = 28 and L = 30) 
a = 0.49(1) or a = 0.51(1) (see Table 2). This is in striking numerical agreement with 
Gribov's formula D^'{k) = {\k\/2)[{k^f + M^]-i/2, which gives a = 0.5. 

3. We have obtained an excellent 2-pole fit for D^'^{k). Our fit indicates that the 
poles occur at complex = x ± iy. The real part x is quite small and compatible with 
0. Remarkably, a pole in at purely imaginary = ± zy agrees with the Gribov 
propagator D^^ik) = {\k\/2)[{Pf + M^]-^/^, with = m^. Note that only a purely 
imaginary pair of poles gives a correction to the free equal-time propagator D^^{k) ~ 
|fc|~^(l — I'^y) °f relative order (A;^)^ with coefficient of dimension (mass)'*. It may not 
be a coincidence that this is the dimension of the gluon condensate (-F^), which is the lowest 
dimensional condensate in QCD. Because of the gauge invariance of the location of the 
poles, by virtue of the Nielsen identities [1], and because of the theoretical suggest ivcness 
of our result, we are encouraged to report the values = ± iy, for y = 0.671(7) in 
lattice units, or y = 0.590(7) GeV^, M = ^ 0.768(4) GeV for the location of the 
gluon poles in k^. 

4. The Coulomb gauge offers a definition of the running coupling constant, ^^coui(^) ~ 
k'^D44{k), which has advantages for numerical determination. It is less subject to fluctua- 
tion than the determination of from the 3-point function in the Landau gauge [25] , but 
the extrapolation in /? remains to be done. 

5. Our data for D44(/c) require a cut at low momentum to eliminate lattice artifacts. 
After this cut, the running coupling constant defined by n^coui(^/^QC7D) = P-D44(fc) 

A similar result, i.e. a transverse gluon propagator that at zero momentum k extrapolates 
to in the limit of infinite lattice volume V, has been recently obtained for pure SU (2) lattice 
gauge theory in the three-dimensional case and in the magnetic sector at finite temperature [26] . 
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extrapolates to low momentum in accordance with infrared slavery, namely the running 
coupling constant 9coui(^/ -^QCd) diverges in the zero-momentum limit. 

6. The observed strong enhancement of the instantaneous color-Coulomb potential 
V{k) and the strong suppression of the equal-time would-be physical gluon propagator 
D^^{k) both at low k, strongly support the confinement scenario of Gribov [6], [7]. In 
addition to this qualitative agreement, we note excellent numerical agreement of our fit 
to Gribov's formula D^'{k) = {\k\/2}[{P)^ +M'^]-^/^, reported in 1, 2, and 3 above. If 
this excellent fit is maintained at larger (3 values, then it appears that we have obtained a 
quantitative understanding of D^^{k). 
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Appendix A. Numerical gauge fixing 

The minimal lattice Coulomb gauge is defined by two gauge-fixing steps. In the first 
step the spatial link variables U^u SU{N), for i = 1, 2, 3, are made as close to unity as 
possible by minimizing the "horizontal" minimizing function 

i^hor,f/(^) = Re TV(1 - ^Usti) , (A.l) 

with respect to gauge transformations g^f The sum extends over all horizontal or space- 
like links, and the minimization is done independently on each time-slice t — X4. (The 
gauge transform ^U^^ of the link variable U^^ is defined by ^Uj^y = U^y Qy, where 
y = X + fi and jl is the unit vector in the fi direction.) The stochastic over-relaxation 
algorithm does not necessarily yield the absolute minimum of the minimizing function but 
leads in general to one of several local minima. Different minima correspond to different 
Gribov copies. For the lattice volumes 16^ and 20"^ we have checked that the dependence 
of the gluon propagators D^^{k) and D44(A;) on which Gribov copy one ends up is of the 
order of magnitude of the numerical accuracy, in agreement with Ref. [27]. 

After this step the lattice gluon field is 3-dimensionally transverse [see Eq. (A. 4) 
below]. But the gauge fixing is as yet incomplete because it leaves a t-dependent but 
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^-independent gauge transformation arbitrary. This arbitrariness is fixed in the second 
step in which the time-like hnk variables Ugt4: ^ire made as cfose to unity as possible by 
minimizing the "vertical" minimizing function 

x,t 

with respect to :r-independent gauge transformations gt- The sum extends over all vertical 
or time-like links. This gauge fixing of the vertical links does not alter the spatial correlator 
D^^ik) nor the color-Coulomb potential V(^*^). However it reduces the non-instantaneous 
part of D44(^), so that it is suppressed compared to the instantaneous part V{k) and 
vanishes with the lattice spacing in the continuum limit. 

This gauge fixing, in which first F^or is minimized and then F^er is minimized is 
equivalent to the limit of the interpolating gauge in which the single function, depending 
on a real positive parameter A, 

F = Fhor + AF,er (A.3) 

is minimized, and the limit A — > is taken. 

With this gauge fixing, the link variables should approach unity in the continuum 
limit in the sense that lim^^oo(l/2)Tr(t/a;^) = 1. In our study at /5 = 2.2 we ob- 
tain {{1 /2)Tr (Uxi)) = 0.86249(3) for the space-like links, whereas for time-like links 
{{l/2)Ti{Ux4)) = 0.221(6) on lattice volume 14^. (The dependence on the volume is 
not strong.) Thus we are quite far from the continuum limit for quantities such as V{k) 
that depend on the vertical links, and they may exhibit significant /^-dependence in the 
extrapolation to the continuum limit. This will be reported subsequently. 

The gauge fixing just described produces a configuration U which is a local minimum 
of Fhor,uig) and Fver,i7(^) at ^ = 1. At a local minimum the minimizing functions are 1) 
stationary under infinitesimal variations SF = and 2) the matrix of second variations 
of the minimizing function is non-negative, S^F > 0. We now comment on implications 
of these properties for gauge-fixed configurations U. At a local minimum, the horizontal 
minimizing function -Fhor,t/(fi') is stationary with respect to infinitesimal variations Qx 
gri;{l+uJx)- Here = f^^x element of the Lie algebra of the SU{N) group, with anti- 
hermitian basis satisfying [t",/:''] = f°'^H'^ and T\:{t°'t^) = —^5°'^. The corresponding 
variation of Fhor is given by 

5F^orA9) = E^[K-+M - ^^t)('UxU - ^4^)] • (A.4) 
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For a configuration U which is a local minimum (at = 1), this quantity must vanish for 
all ujr^tj which gives 

J2iMt,i - Ax-i,t,i) = • (A.5) 

i 

Here is the lattice gluon field defined by 

Al^^-Tr[nU,,-Ul)l (A.6) 

which is a lattice analog of the continuum connection Equation (A.5) is the lattice 

transversality condition for spatial directions, which is the defining condition for the lattice 
Coulomb gauge. 

Since the gauge-fixed configuration f/ is a local minimum of Fhor,J7(5') (at gx = 1), its 
second variation is non-negative 

S^Fhor,u{g) = {ou, M{U)u) > for aU u; . (A.7) 

Here M(U) is the lattice Faddeev- Popov matrix defined on a given time-slice t by 

(w,M(?7)a;) = -^^I¥[(w-^. -a;af)(?7a-ic^5r+i-c^xt^a-i^ , (^^.8) 

and we have suppressed the index t which is common to all variables. The positivity of 
M{U) is a condition on configurations U which, together with the transversality condition 
(A.5), defines the lattice Gribov region, whose boundary is the Gribov horizon. 

[That this is a highly restrictive condition is suggested by the following considera- 
tion. The Faddeev-Popov matrix M{U) for SU{N) gauge theory is a symmetric matrix 
of dimension V{N'^ — 1), where V is the (large) number of sites of the lattice, so it has 
V{N'^ — 1) eigenvalues. Configuration space is divided into V{N'^ — 1) -|- 1 different regions 
Rn according to the number n = 0, . . . , V{N'^ — 1) of positive eigenvalues of M{U). Of 
these, the Gribov region consists of the single region i2y(-jv2_i) that includes U^i = 1. For 
the SU{2) group at least, all regions are populated. To see this, observe that for Uxfj, = 1, 
we have M(l) = —A, whereas for = ~1, we have M(— 1) = A, where A is the lattice 
Laplacian. In these 2 cases, depending on the sign, the configuration U = ±1 is in region 
Rv(N'2-i) or Rq. By continuity therefore all V{N'^ — 1) + 1 different regions are populated. 
Similar considerations apply to Fver,c/(5')-] 
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Table Captions 



Table 1. Fit of the transverse gluon propagator D^^{k) using Eq. (2.8). The resulting 
fitting parameters, x^/d.o.f. and goodness-of-fit Q are reported for each lattice 
side L. 

Table 2. Extrapolation to infinite lattice side L for the five fitting parameters appearing in 
Eq. (2.8) (we always have r = 0). Five difi^erent cases are considered (see Section 
2). 

Table 3. Fit of the running coupling constant ^coui(^) = k^D44{k) using Eq. (2.11). The 
resulting fitting parameters, x^/d.o.f. and goodness-of-fit Q are reported for three 
different values of the low-momentum cut kl,:„. 
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L 


z 


r 


a 


X 




xVd.o.f. 


Q 


14 


3.93(33) 


3.86(35) 


0.76(5) 


0.36(6) 


1.07(5) 


0.42 


73.9% 


16 


4.79(42) 


2.72(28) 


0.70(5) 


0.27(8) 


1.08(7) 


1.46 


21.1% 


18 


5.91(53) 


2.27(33) 


0.57(7) 


0.38(13) 


1.11(13) 


3.45 


0.4% 


20 


5.78(19) 


1.80(8) 


0.66(2) 


0.19(4) 


0.99(4) 


0.45 


84.5% 


22 


5.81(23) 


1.57(8) 


0.67(3) 


0.09(4) 


0.97(5) 


0.59 


76.5% 


24 


6.33(25) 


1.43(8) 


0.62(3) 


0.12(4) 


0.99(6) 


0.86 


55.0% 


26 


6.99(22) 


1.08(5) 


0.62(2) 


0.12(3) 


0.82(4) 


0.60 


79.8% 


28 


7.21(28) 


1.05(6) 


0.60(3) 


0.14(4) 


0.82(5) 


0.72 


70.6% 


30 


7.08(45) 


0.93(8) 


0.64(4) 


0.06(6) 


0.78(8) 


1.54 


11.0% 



Table 1 
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p 

min 


# data points 


z 


a 


X 




xVd.o.f. 


Q 






1 7 OfOdQ) 












0.5 


23 


12.8(20) 


0.45(7) 


0.11(17) 


0.45(6) 


0.98 


48.1% 


0.3 


25 


13.4(11) 


0.43(4) 


0.17(8) 


0.42(3) 


0.94 


53.8% 


0.5 


23 


11.5(2) 


0.49(1) 





0.47(4) 


0.96 


50.9% 


0.3 


25 


11.1(1) 


0.51(1) 





0.39(2) 


1.22 


21.7% 



Table 2 



^min # data points ^loui ^ ^' = Ti^ x^/d-o.f. Q 



0.3 25 0.84(8) 1.14(5) 0.165(9) 2.2(1) 0.59 92.8% 

0.5 23 1.0(2) 1.2(1) 0.18(2) 1.9(3) 0.44 98.3% 

1.0 20 1.2(9) 1.3(4) 0.19(7) 1.7(7) 0.48 95.8% 



Table 3 
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Figure Captions 



Fig. 1. Plot of the gluon propagators D^^{k) and -D44(A;) as a function of the square of 
the lattice momentum /c^ for L = 28 (symbols * and Q respectively) and L = 30 
(symbols A and V respectively). Notice the logarithmic scale in the y axis. Error 
bars are one standard deviation. 

Fig. 2. Plot of the gluon propagator D^^{k) as a function of the square of the lattice 
momentum P for L = 14(*), 20(A) and 26(0) fits of these data using Eq. 
(2.8) with the parameter values reported in Table 1. Error bars are one standard 
deviation. 

Fig. 3. Plot of the gluon propagator D^'^(k) as a function of the square of the lattice 
momentum P for L = 16(*), 22(A) and 28(0) fits of these data using Eq. 
(2.8) with the parameter values reported in Table 1. Error bars are one standard 
deviation. 

Fig. 4. Plot of the gluon propagator D^'^(k) as a function of the square of the lattice 
momentum k'^ for L = 18(*), 24(A) and 30(0) fits of these data using Eq. 
(2.8) with the parameter values reported in Table 1. Error bars are one standard 
deviation. 

Fig. 5. Fit of the parameter z{L) (see Table 1) as a function of 1/L using a — b/L^. 
We obtain a = 17.2(249), b = 35.9(117), c = 0.38(83), xV^-o.f. = 0.89 and 
goodness-of-fit Q = 50.1%. 

Fig. 6. Fit of the parameter a{L) (see Table 1) as a function of 1/L using a+log(l + 6/L). 

We obtain a = 0.48(5), b = 4.0(13), xV^-O.f. = 0.66 and goodness-of-fit Q = 
70.6%. 

Fig. 7. Fit of the parameter x{L) (see Table 1) as a function of 1/L using a + log(l + 6/L). 

We obtain a = -0.2(1), b = 10.0(34), xV^-O-f. = 0.94 and goodness-of-fit Q = 
47.4%. 

Fig. 8. Fit of the parameter y'^{L) (see Table 1) as a function of 1/L using a+log(l+6/L). 

We obtain a = 0.45(1), b = 13.6(45), xVd-O-f- = 1-05 and goodness-of-fit Q = 
39.3%. 

Fig. 9. Fit of the product {rz){L) (see Table 1) as a function of 1/L using exp(a)/L^. We 
obtain a = 5.6(4), b = 1.1(1), xV^-o.f. = 0.36 and goodness-of-fit Q = 92.6%. 

Fig. 10. Fit of the zero- momentum transverse gluon propagator D*^(0, L) as a function of 
1/L using exp(a)/L''. Considering the data for L > 16, we obtain a = 2.43(4), 
b = 0.50(1), xVd-o.f. = 0.69 and goodness-of-fit Q = 65.6%. 
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Fig. 11. Plot of the gluon propagator D^'^{k) as a function of the square of the lattice 
momentum k'^ for L = 28(*) and 30(A) and fits of these data using Eq. (2.8) 
with the parameter values and low momentum cut-offs reported in the fourth 
(solid line) and fifth (dotted line) rows of Table 2. Error bars are one standard 
deviation. 

Fig. 12. Plot of the running coupling constant xf^coui(^) = ^^-044 (^) as a function of the 
square of the lattice momentum k'^ for L = 28(*) and 30(A) and fits of these data 
using Eq. (2.11) with the parameter values and low momentum cut-offs reported 
in the second row of Table 3. Error bars are one standard deviation. 
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